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Abstract 

We study strict local martingales via /i-transforms, a method which 
first appeared in Delbaen-Schachermayer. We show that strict local 
martingales arise whenever there is a consistent family of change of 
measures where the two measures are not equivalent to one another. 
Several old and new strict local martingales are identified. We treat 
examples of diffusions with various boundary behavior, size-bias sam- 
pling of diffusion paths, and non-colliding diffusions. A multidimen- 
sional generalization to conformal strict local martingales is achieved 
through Kelvin transform. As curious examples of non-standard be- 
havior, we show by various examples that strict local martingales do 
not behave uniformly when the function (x — K)'^ is applied to them. 
Implications to the recent literature on financial bubbles are discussed. 

1 Introduction 

Local martingales which are not martingales (known as "strict" local martin- 
gales) arise naturally in the Doob-Meyer decomposition and in multiplicative 
functional decompositions, as well as in stochastic integration theory. They 
are nevertheless often considered to be anomalies, processes that need to be 
maneuvered by localization. Hence, studies focussed purely on strict local 
martingales are rare. One notable exception is the article by Elworthy, Li, & 
Yor [9] who study their properties in depth. On the other hand, applications 
of strict local martingales are common. See the articles by Bentata & Yor 
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[2], Biane & Yor [3j, Cox & Hobson [3], Fernholz & Karatzas [12j, and the 
very recent book-length preprint of Profeta, Roynette, and Yor [33j. 

Our goal in this paper is to demonstrate that strict local martingales 
capture a fundamental probabilistic phenomenon. Suppose there is a pair 
of probability measures on a filtered probability space, where one strictly 
dominates the other in the sense that their null sets are not the same. Then 
one can establish a correspondence between pairs of continuous processes 
such that one is a martingale under one measure and the other is a strict 
local martingale under the changed measure. This has been made precise in 
Propositions [T] and [6j The continuity of the time parameter is crucial since 
there are no strict local martingales in discrete time. 

The above idea is not new and has appeared in the literature under 
various disguises. For example, for positive local martingales, such a phe- 
nomenon was identified by Delbaen and Schachermayer in [6]. We, however, 
do a systematic exploration to identify several strict local martingales as an 
application of the previous result, in diverse topics such as diffusions con- 
ditioned to exit through a subset of the boundary of a domain, size-biased 
sampling of diffusion paths, and non-colliding diffusions such as Dyson's 
Brownian motion from the Random Matrix Theory. In addition, we show 
how usual stochastic calculus with respect to strict local martingales can be 
made more friendly by incorporating the change of measure. Our results in 
this direction are related to recent work by Madan and Yor |27j . 

We also prove a multidimensional analogue of our results where a strict 
local martingale in one-dimension is replaced by a conformal local martingale 
in three or more dimensions where at least one coordinate process is strict. 
The analysis exploits the Kelvin transform from classical potential theory. 

The name h-transforms in the title refers to the classical technique (orig- 
inally due to Doob), where the law of a Markov process is changed by using 
a nonnegative harmonic function. Aided by Girsanov's theorem, this usu- 
ally leads to a description of the process conditioned to stay away from the 
zero-set of the harmonic function. As our examples will demonstrate, this 
is a fertile set-up for a plethora of strict local martingales. For more details 
about classical /i-transforms, please see Revuz and Yor p. 357]. 

The rest of the paper analyzes the effect of convex functions applied to 
strict local martingales. Although these are always local submartingales, 
the expectations of the processes can be of arbitrary shapes. Similar exam- 
ples were also identified in [9]. In particular, this has some consequences 
in mathematical finance. A natural source for local martingales in math- 
ematical finance is the condition of No Free Lunch with Vanishing Risk 
(see Delbaen and Schachermayer [5]). Roughly, it states that in a financial 
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market the no arbitrage condition is equivalent (in the case of continuous 
paths) to the existence of an (equivalent) "risk neutral" probability measure 
Q which turns the price process into either a martingale or a strict local 
martingale. The implications of our results can be readily understood if we 
assume that the risk neutral measure produces a (one-dimensional) price 
process X = {Xt)t>o that is a strict local martingale. In that case, for ex- 
ample, the process Yt = {Xt — K)^ need not be a submartingale, and the 
function 1 1— >• E{{Xt — K)^} need no longer be increasing, contradicting the 
usual wisdom in the theory. The original purpose of this paper was to un- 
derstand this phenomenon better, motivated in particular by the role local 
martingales play in financial bubbles (cf [20] and [21j). 

2 A method to generate strict local martingales 

Let ^il, {J-t}j>Q^ be a filtered sample space on which two probability mea- 
sures P and Q are defined. We assume that P is locally strictly dominated 
by Q, in the sense that P is absolutely continuous with respect to Q (-P <C Q) 
on every (see [Ml Chap. VIII]). Let ht denote the Radon-Nikodym process 
dP/dQ on the a-algebra J^f Assume that 



where tq = inf {s > 0, hs = 0}. We claim the following result. 

Proposition 1. Assume that h is a continuous process Q-almost surely. 
Let {ft, t > 0} be a continuous Q-martingale adapted to the filtration {Ft}- 
Suppose either of the two conditions hold: 

(i) f is uniformly integrable, E^^fo) ^ 0, and Q (tq < oo) = 1, or, 
(a) f is nonnegative, Q{ao > tq) > 0, where ao = mf{s > : fs = 0}. 

Then the process Nt := ft/ht is a strict local martingale under P. However, 
if f is nonnegative and Q{cfq > tq) = 0, the process Nt is a martingale. 

Proof of Proposition [H Let us first show that is a local martingale. Con- 
sider the sequence of stopping times 



Then, it follows that PiliiUk^ooTk = oo) = 1 (see [Ml P- 326]). Thus, for 
any bounded stopping time r, we get by the change of measure formula 



Q (to < oo) > 0, 



(1) 



Tk:=mf{s>0: hs<l/k}, k=l,2,.... 
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The final equality is due to the fact that / is a Q-martingale and r ATk is 
a bounded stopping time. 

Since NrArk has the same expectation for all bounded stopping times r, 
it follows that A^-at^. is a martingale. The local martingale property now 
follows. 

To show now that it is not a martingale, we compute the expectation of 
A^^. Note that h is a nonnegative Q-martingale. Thus, again applying the 
change of measure formula, we get 

[Nt) = eQ {ftl{r,,>t}) , (2) 

where tq is the hitting time of zero for h. 

Now suppose (i) {ft} is a uniformly integrable martingale and Q (tq < oo) = 
1. By uniformly integr ability we see from ^ that 

lim E^{Nt) = 

i— >oo 

which shows that N is not a martingale, since E^{Nq) = E^{fQ) is assumed 
to be non-zero. 

Finally suppose (ii) Q{(To > tq) > holds. Since zero is an absorbing 
state for the nonnegative martingale /, by ([T]) there is a time t > such 
that Q{{ft > 0} n {tq < t}) > 0. For that particular t, we get from ^ that 

E"" (Nt) < E'^ift) = E^ifo) = E^'iNo). 

This again proves that the expectation of Nt is not a constant. Hence it 
cannot be a martingale. 

The final assertion follows from ([2]) by noting that A^, by virtue of being 
a nonnegative local martingale, is a supermartingale which has a constant 
expectation. Thus it must be a martingale. □ 

It is clear that some condition is necessary for the last theorem to hold 
as can be seen by taking ft = ht which results in a true martingale Nt = 1. 

Several examples of old and new strict local martingales follow from the 
previous result by suitably choosing a change of measure and the process /. 
We describe some classes of examples below. 

2.1 Diffusions with different boundary behaviors 

One of the earliest uses of change of measures was to condition a diffusion 
to exhibit a particular boundary behavior. This typically involves a change 
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of measure that is locally strictly dominated. We give below two examples 
to show how local martingales arise from such a set-up. 

Example 1. The 3-dimensional Bessel process, BES(3), (see Karatzas & 
Shreve [231 page 158] for the details) is the (strong) solution of the stochastic 
differential equation: 

dXt = ^dt + d^t, Xo = xo> 0, (3) 

where f3 is an one-dimensional standard Brownian motion. It is also the 
law of the Euclidean norm of a three-dimensional Brownian motion. The 
reciprocal of this process, known as the inverse Bessel process, serves as a 
prototypical example of a local martingale which is not a martingale. 

This is an immediate example of Proposition [1] where we consider the 
canonical sample space for a Brownian motion, define Q to be law of the 
Brownian motion, starting at xq, and killed at zero, and let ht be the coor- 
dinate map LOt/xQ. Finally we take ft = 1- Details on why Q is the law of 
the Bessel process are well-known and can be found in [34, p. 252]. 

Several other examples of the same spirit can be derived. In particular, 
for any Bessel process X of dimension 6 > 2, it is known that X^~^ is a 
strict local martingale. This can be proved similarly as above using the 
martingale y^~^, where y is a Bessel process of dimension (4 — 6) (which 
can be negative) absorbed at zero. See the article ^ for the details. 

Example 2. The previous example can be easily extended to a multi- 
dimensional form. Let D be an open, connected, bounded domain in M" 
(n > 2) where every point on the boundary is regular (in the sense of \23\ 
p. 245]). Consider an n-dimensional Brownian motion X starting from a 
point xq G D getting absorbed upon hitting the boundary of D, say dD. 
Let Q denote the law of the process {XtAm, ^ — ^i- Since D is bounded, 
Q{td < oo) = 1. For any bounded measurable nonnegative function u on 
dD, one can construct the following function 

f{xo) = E'^{u{Xr^)\Xo = xo). (4) 

In fact, / is the solution of the Dirichlet problem on D with boundary data 
u. By the Markov property, it follows that f{Xt/\Tjy) is a martingale. 

Let Bi be a connected (nontrivial) proper subset of the boundary dD, 
and let the function u be one on Bi and zero elsewhere, i.e, u{x) = If^^^Bi)- 
In that case, the resulting harmonic function in D is given by 

v{xo) = Q {Xrjj £ Bi\ Xo = Xq) . 
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Then v{xo) > 0, and by the tower property, h{Xt/\Ta) = v{Xt/\Tjj) /v{xo) is 
a nonnegative martingale starting from one. Let P denote the change of 
measure using the Radon-Nikodym derivative process h^Xt^To)- This gives 
us a process law which can be interpreted as Brownian motion, starting from 
xq, conditioned to exit through Bi. 

As a corollary of Proposition [H we get the following result. 

Proposition 2. Let u he any nonnegative function on dD such that u > on 
some subset of dD\Bi of positive Lebesgue measure. Let f be the harmonic 
extension of u given by ([4]) . Then the process 

* viXt/,rr,) 
is a strict local martingale under P. 

Proof. Follows from Proposition [T] by verifying condition (ii) since u^X-nj ) 
need not always be zero when X^-^ ^ Bi. □ 

2.2 Size-biased sampling of diffusion paths 

Another class of interesting examples follow in the case of the size-biased 
change of measure. Size-biased sampling has been often discussed in con- 
nection with discrete distributions, see, for example, the article by Pitman 
|30j and the references within. It usually involves a finite or countable col- 
lection of numbers {pi,P2,P3, • • •} such that each pi > and "^^Pi = 1- A 
size-biased sample from this collection refers to a sampling procedure where 
the sample p has the distribution 

P{p = Pi)=Pi, foraUz = l,2, 

One can now remove this chosen sample from the collection, renormalize it, 
and repeat the procedure. This is closely connected to urn schemes where 
each Pi refers to the proportion of balls of a color i that is in an urn. If 
one randomly selects a ball, the color of the chosen ball has the size-biased 
distribution. 

One can similarly develop a concept of size-biased sampling of diffusion 
paths as described below. Consider n non- negative diffusions {^i, . . . ,Xn} 
running in time. Fix a time t, and look at the paths of the diffusions 
during the time interval [0,t]. Denote these random continuous paths by 
Xi [0, t] , X2 [0,t], . . . , Xn [0, t] . Sample one of these random paths with prob- 
ability proportional to the terminal value Xi{t). That is, the sampled path 
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has the law 



Y[0,t] = Xi[0,t\, with probability 



Xi{t) 



How can one describe the law of y ? In general the law of Y might not 
be consistently defined as time varies. Nevertheless there are cases where it 
makes sense. We show below the example when each Xi is a Bessel square 
process of dimension zero (BESQ°). This is the strong solution of the SDE 



where /? is a one-dimensional Brownian motion. BESQ is a nonnegative 
martingale also known as the Feller branching diffusion, in the sense that it 
represents the total surviving population of a critical Galton- Watson branch- 
ing process. Indeed, our treatment here of size-biased transforms of BESQ 
processes is inspired by size-biased transforms of Galton- Watson trees (see 
the work by Lyons- Pemantle- Peres [26] and references to the prior literature 
referred to there). 

The interpretation of such dynamic size-biased sampling when every Xi is 
a BESQ'' is straightforward. When each Xi has the law of Feller's branching 
diffusion, they represent a surviving population from n critical branching 
processes. We can do size-biased sampling at different time points from 
these populations. The construction below describes the joint law of these 
samples. 

Consider the canonical sample space for multidimensional diffusions, i.e., 
the n-dimensional continuous path space C"[0,oo), coupled with the usual 
right-continuous and complete (with respect to the diffusion we describe be- 
low) filtration generated by the coordinate maps. Let ojt = {'^ti'^)-, ■ ■ ■ , ^t{n)) 
denote a sample path. Also, as before, we take ti;[0,t] to denote the path 
during time-interval [0, t]. Let Q be the joint law of n independent BESQ'' 
processes starting from positive points (^(l), . . . , 2;(n)). To keep matters 
simple, we assume all the z{i)' s are the same and equal to z. This induces 
the following exchangeability property. 

Let Tj : C"'[0,oo) — )• C"[0,cxd) be a permutation operator on the sample 
space such that Tjw(l) = uj{i) and 

TiUj{j + 1) = w(j), if j < i, and, Tiw(j) = uj{j), if j > i. (5) 

Thus, Ti puts the ith coordinate as the first, and shifts the others appropri- 
ately. By our assumed exchangeability under Q, each TiU has also the same 




z > 
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law Q. To define the size-biased change of measure we need the following 
lemma. 



Lemma 3. For any i, the process 



Mt{i) 



is a martingale under Q. 



ujt{l) + ... +ujt{n)'- 



t > 0, 



Proof. Without loss of generality take i = 1. Let Z{1), Z{2), . . . ,Z{n) be 
independent BESQ°. The sum C = Z{1) + ... + Z{n) is another BESQ° 
process. Thus by Ito's rule we get 
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d{Ztil)/Ct) = Zt{l)d(- ]+-dZt{l)+d(Ztil),Q 



CtJ Ct 

local martingale +^4C.dt-i^d. 



Ct 



This proves that the ratio process is a local martingale. But since it is 
bounded, it must be a martingale. □ 

Define the size-biased sampling law P on (C"[0, oo), {Tt}) by 



P{A) = 



^^u}t{l) + ---+ujt{n) 



{TiUj[0,t\&A} 



, ioxaWAeFt. (6) 



Note that, by Lemma [3l P defines a consistent probability measure on the 
filtration {Ft}. We show below that P is strictly locally dominated by Q 
and compute the Radon-Nikodym derivative. 

Recall that, under Q, each Tjo; has the same law Q. Thus, we can 
simplify expression ([6]) to write 



P{A) = nE^ 



-1 



{u;[0,i]GA} 



for all Ae Ft. 



H h ujt{n) 

This proves that P Q and the Radon-Nikodym derivative is given by 

nutil) 



ht 



ujt{l) H \-uJt{n)' 



Since under the BESQ'' law, every coordinate can hit zero and get absorbed, 
the above relation is a locally strict domination and hence leads to examples 
of strict local martingales. 
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It can be verified, using Girsanov Theorem, that under P, the coordinate 
processes remain independent. The first coordinate w(l) has law BESQ 
of dimension 4, and the rest of the coordinates are still BESQ*^. These 
correspond to the spine decomposition of critical GW trees, conditioned on 
surviving forever, as in [26j . 

As an immediate corollary of Proposition [1] we get the following result. 

Proposition 4. Let {Z{1), . . . ,Z{n)) be continuous processes whose law is 
the size-biased sampled BESCf law P described above. Let Q = Zt{l) + . . . + 
Zt{n) be the total sum process. Then the processes 

are all strict local martingales. However, the process 

n 

Mt = CtYlZt{i), t>0, 

1=2 

is a true martingale. 

Proof. Let Z(l), . . . , Z{n) be iid BESQ^ processes starting from z > 0. Then 
C, Z{2), and YYi=2^ti^) true martingales which can remain positive 

when Zt{l) = 0. The result now follows from Proposition [U condition (ii). 

On the other had, the process Yli=i -^4(0 is a martingale that always 
hits zero before Zt{l). Thus, M is a true martingale. □ 



2.3 Non-colliding diffusions 

Our third class of examples are cases when the change of measure leads 
to non-intersecting paths of several linear diffusions. Probably the most 
important example of this class is Dyson's Brownian motion which is a 
solution of the following n-dimensional SDE: 

dXt{i) = y] ■ ^ , dt + dBt{i), t>0, i = l,...,n. (7) 

Here (-B(l), • • • ,B{n)) is an n-dimensional Brownian motion. It appears in 
the context of Random Matrix Theory. Please see the survey article by 
Konig [21] and the original paper by Dyson [7] for the details (including the 
definition of the Gaussian Unitary Ensemble) and the proofs. 
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Theorem 1. For any i = 1,2, ... ,n, and i < j < n, let {Mt{i, i), t > 0}, 
{M^{i,j), t > 0}, and {M^{i,j), t > 0} be independent real standard 
Brownian motions, starting at zero. The Hermitian random matrix Mt = 
{Mt{i,j), l<i,j <n), with 

Mt{i,j) = M^{i,j) + iMl{i,j), i < j, 

has the distribution of the Gaussian Unitary Ensemble at time t = 1. Then 
the process (At, t >0) ofn eigenvalues of Mt satisfies SDE ([7|). It can inter- 
preted as a conditional Brownian motion in M", starting at zero, conditioned 
to have 

At(l) < At(2) < • • • < \t{n), for all t > 0. 

What is interesting is that the process in ([7]) can be obtained as a strict 
local domination relation from the n-dimensional Wiener measure using the 
harmonic function 

^n{x)= W {x{j)-x{i)), x = {x{l),...,x{n)), 

l<i<j<n 

which is the well-known Vandermonde determinant. That the function A„ 
is harmonic can be found in [24, p. 433] where it is shown that A„(Wj) is 
a martingale when W is an n-dimensional Brownian motion and that the 
law of the process in ([7| (say P) can be obtained from the n-dimensional 
Wiener measure Q, by using ht = An{Wt) as the Radon-Nikodym derivative. 
Since An{W) is zero whenever any two Brownian coordinates W{i),W{j) 
are equal ("collide"), we are in the scenario of Proposition [TJ We prove the 
following result. 

Proposition 5. Consider Dyson's Brownian motion (At(l), . . . , At(n)) in 
d?!, and for some m < n consider the process 

l<i<jr<m 

Then the process Nt = Ara{t)/An{t) is a strict local martingale. As a con- 
sequence if we consider the Vandermonde matrix-valued process 

At = iXl-\{}, i,j = 1,2,..., n), t>0. 

Then every process A^^{n, i) is a strict local martingale for i = 1,2, . . . ,n. 
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Proof. Note that, if W is an n-dimensional Brownian motion then Am{Wt), 
for any m < n, is a true martingale. Since m < n, it is possible to have 
AfniWt) to be positive when An{Wt) = 0. The result that is a strict local 
martingale now follows from Proposition [TJ 

For the second assertion we use adjugate (or, classical adjoint) relation 
for the inverse. Let det(i?) for a square matrix B refer to its determinant. 
Then 

A^Hn, i) = ^et(^) l^*^*' ' ^ = 1' 2, . . . , n, 

where At{i,n) is the matrix obtained from At by removing the ith row and 
the nth column. 

Now, At is the Vandermonde matrix, so its determinant is equal to A„(f). 
If we remove the ith. row and the nth column from At we get an (n — 1) x (n — 
1) order Vandermonde matrix of all the Aj's except the ith. Its determinant 
is again a Vandermonde determinant. Each of the ratios 

det (it(i,n)) 
det{At) ' 

is a strict local martingale by our earlier argument (for m = n — 1). This 
completes the proof. □ 

As a final area of applications of /i-transforms which are locally strictly 
dominated, let us mention the theory of measure-valued processes, in par- 
ticular, the Dawson- Watanabe superprocesses. A well-known example is 
conditioning a superprocess to survive forever, which can be done by chang- 
ing the law of a superprocess by using the total mass process (which is a 
martingale) as the Radon-Nikodym derivative. Please see the seminal article 
in this direction by Evans & Perkins [iT], the book by Etheridge [10], and a 
follow-up article on similar other /i-transforms by Overbeck [29]. 

3 A converse to the previous result 

As a converse to Proposition [1] it turns out that all strict local martingales 
which remain strictly positive throughout can be obtained as the reciprocal 
of a martingale under an /i-transform. This was essentially proved by Del- 
baen and Schachermayer [6] in 1995 in their analysis of arbitrage possibilities 
in Bessel processes. We replicate their theorem below. The construction is 
related to the FoUmer measure of a positive supermartingale |14j . 
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Before we state the result we need a technique which adds an extra ab- 
sorbing point, oo, to the state space originally inspired by the work of 
P. A. Meyer [2H]. We follow closely the notation used in [6]. The space of 
trajectories is the space Coo[0, T] or Coo[0, oo) of continuous paths to defined 
on the time interval [0, T] or [0, oo) with values in [0, oo] with the extra prop- 
erty that if uj{t) = oo, then uj{s) = oo for all s > t. The topology endowed is 
the one associated with local uniform convergence. The coordinate process 
is denoted by X, i.e., X{t) = uj{t). 

Theorem 2 (Delbaen and Schachermayer, Theorem 4 in ^ ). If R is a 

measure on C[0, 1] such that X is a strictly positive strict local martingale, 
then 

(i) there is a probability measure R* on Cod[0,1] such that M = 1/X is 
an R* martingale. 

(a) We may choose R* in such a way that the measure R is absolutely 
continuous with respect to R* and its Radon-Nikodym derivative is 
given by dR = MidR* . 

The following result is a corollary. 

Proposition 6. Let R be a probability measure on C[0, oo) under which the 
coordinate process X is a positive strict local martingale starting from one. 
Then there exists a probability measure Q on the canonical space such that 
X is a nonnegative martingale under Q and the following holds: 

(i) The probability measure defined by 

P{A) := E"^ {XdA) , V^GJi, t>0, (8) 
is the law of the process {1/Xt^ i > 0} under R. 

(ii) X is a strict local martingale if and only if Q{tq < oo) > 0, where 
tq = inf {t > : Xf = 0} is the first hitting time of zero. 

Proof. Let us first construct a probability measure Qi on C[0, 1] such that 
([8]) holds for < t < 1. To wit, consider the restriction of R on C[0, 1]. The 
coordinate process, under the restricted measure, remains a positive strict 
local martingale. Thus, there is an R* such that Theorem [2] holds. Define 
Qi to be law of the process {Mt, < t < 1} under R*. Then, we claim 
that the finite dimensional distributions of the coordinate process X under 
Qi is the same as the law of 1/X, under R. To see this, consider m G N 
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many points < ti < t2 < ■ ■ ■ tm-i < < 1- For any choice of nonnegative 
measurable function g : M™ — )• M"*", we get 

E« [g {X^;\X,-\..., XiJ] = E""* [Ml 5 (Mt, ,Mt,,...,Mt^ )] 

= i?«^ [X^g{Xt,,Xt,,...,XtJ]. 

The first equahty above is due to part (ii) of Theorem [2l This proves ^ for 
<t < 1. Since a measure on C[0, 1] is determined by its finite dimensional 
distributions, we are done. 

For any other time t, consider the scaled process {Xm, < u < 1} 
under R. It remains a strictly positive strict local martingale, now running 
between time zero and one. As in the previous paragraph one can construct 
a measure, such that, for any choice of < ti < t2 < ■ ■ ■ < tm < t, we 
have 

E^[g{Xi^\Xi^\...,Xi^)]=E^'^ [X^g{Xt,/t,X,^/„...,Xt^/t)]. (10) 

Let Qt denote the law on C[0, t], such that the law of the scaled coordinate 
process {Xfu, < n < 1} under Qt is fi\. Once we demonstrate that this 
tower of probability measures {Qt, t > 0} is consistent, it follows from 
standard arguments that they induce a probability measure Q on the entire 
space C[0, oo) such that Q restricted C[0,t] is Qt- 

To see consistency, let s < n be two positive numbers, and we show 
that Qu restricted to Ts is Qs- As before, choose m G N, and consider a 
nonnegative measurable function g : M*" — t- M"*". Choose 

< ti < t2 < ■ ■ ■ tm-l < tm < S < U. 

Since X is a martingale under Q^ (and Qs), we have 

eQ- [Xt^g {Xt„..., XtJ] = E'i- [X^g (Xt,,..., XtJ] = 
E^'i [X^g {Xt,,u, Xt^,^)] = E^ [g {X,~\ XfJ] , by m, 

= [Xig {Xt^/,, . . . , Xt^/s)] , again by m, 

= E^^ [Xsg {Xt„..., Xt^ )] = E^' [Xt^g {Xt„ . . . , XtJ] . 

Since Xt^g (^ti, • • • , XtJ is another function of finite dimensional distribu- 
tions, the above identity shows the following. For any function h : (M"*")™" — t- 
M+ such that 

d 

h(x) = 0, if Xm = 0, and < 7: — h(x) < 00, (11) 

oxm 0+ 
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one has 

[h {Xt„...,Xtj] = [h {Xt„...,Xtj]. 

However, under both Q" and Q'^, the coordinate process is nonnegative 
almost surely. For nonnegative random variables the class of functions in 
(jlip is a distribution determining class since indicators of open sets in (M^)™ 
can be approximated by smooth functions satisfying (|lip . This shows con- 
sistency and part (i) of the Proposition is proved. 

For part (ii) we use the change of measure in ^ to claim that 

E^Xt = E^l/Xt = E'^ (xtl{Xt > 0}^^ = Q (Xt > 0) = Q (tq > t) . 

If Q{tq < oo) > 0, then E^Xt is not constant over t which prevents it form 
being a martingale. This completes the proof of the Proposition. □ 



4 Stochastic calculus with strict local martingales 

The advantage of Proposition [6] is that it allows us to transport stochastic 
calculus with respect to strict local martingales to that with actual martin- 
gales via a change of measure. 

Proposition 7. Let R and Q he two probability measures on C[0, oo) such 
that the coordinate process X is a positive strict local martingale under R, 
is a nonnegative martingale under Q, and the relationship ([8]) holds. 

Let Tq denote the hitting time of zero of the coordinate process Xt . Con- 
sider any nonnegative function h : (0, co) — )■ M+. For any hounded stopping 
time T, we get 

E''{h{Xr)) = E'^g{Xr)l{r,>ry (12) 

Here g is the function g{x) = xh{l/x), for all x > 0. 

Now suppose lim-j;_!.o 5(2;) = rj < 00. Define a map g : [0,oo) M by 
extending g continuously, i.e., g{x) = g{x) for x > 0, and g{0) = rj. Then 
we have 

{h{Xr)) = E^giXr) - 7]Q {to < r) . (13) 
Proof of Proposition [71 From ([8|) one gets 

E^'hiXr) = E'^Xrh (^^) 1{,„>,| = E'^g{Xr)l{ro>r}, 

which proves (fT2]) . 
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For the second assertion note that for any nonnegative path co which 
gets absorbed upon hitting zero, the following is an algebraic identity: 

5(Wr)l{ro>r} = ^K) " f?l{ro<r}- 

In particular this identity holds pathwise when u is a. path of a nonnegative 
martingale. Thus we obtain 

i?«<7(^r)l{.o>r} = E'^giXr) - 7?Q(ro < r). 

This proves (fTHj) . □ 

As an example note that when h{x) = {x — a)^ for some a > 0, we get 

{Xr - a)+ = E^{1- aXr)~^ - Q(to < r). (14) 



We have the following corollary. 

Corollary 1. Let h : (0,oo) — )• (0, oo) be a function which is suhlinear at 
infinity, i.e., 

lim ^ = 0. (15) 
Then, for all hounded stopping times t, one has 

E^h{Xr) = E'^g{Xr), g{x) = xh{l/x), x > 0, ^(0) = 0. (16) 
Proof. The first part follows directly from ()13p since 

Tj = lim xh \ — \ = lim — 



\X J x-^oo X 

The second conclusion is obvious. □ 



The previous corollary has some interesting consequences. For example, 
when h is convex, it is not difficult to verify that so is g. And thus both 
h{X) and g{X) are submartingales (under R and Q respectively) by (fT6|) . 
This is in spite of the strictness in the local martingale property of the 
coordinate process under R. Additionally, if h is symmetric with respect to 
inverting x, i.e. h{x) = xh{l/x), then g = h. Hence the strictness of local 
martingales has no effect when these functions are applied. For example, 
E^y/X^ = E^^r- 

Before we end, let us mention that similar results can be obtained from 
the semimartingale decomposition formulas of Madan and Yor [27] and by 
PDF methods, as in Fkstrom and Tysk [8j. 
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5 Convex functions of strict local martingales 



Strict local martingales are known for odd behavior which is not shared by 
martingales. For example, a convex function of a martingale is always a 
submartingale. This need not be the case with local martingales. However 
if A'^ is a nonnegative strict local martingale, and /i is a convex function 
sublinear at infinity, then h{Nt) is again a submartingale. This is in contrast 
to functions which are linear at infinity. For example, in the case of h(x) = x, 
the process is actually a supermartingale. Here we demonstrate another 
example, for the function (x — K)~^ with K > applied to the inverse Bessel 
strict local martingale in Example 1. For a complementary result, see |35| . 



Proposition 8. Let Xt he a BES(3) process starting from one. For any 
real K G [0,1/2], the function t i— )■ E{{l/Xt — K)'^} is strictly decreasing 
for all t £ (0,oo). However, if K > 1/2, the function t H> E{{l/Xt - K)+} 
is increasing in a neighborhood of zero but is strictly decreasing for all large 
enough t. In particular, it is strictly decreasing for 



2K - 1 

Proof of Proposition [HI Let i? be a one-dimensional Brownian motion start- 
ing from one and absorbed at zero. From the change of measure identity in 
we get 

h{t) := E{{l/Xt - K)+} = Eil - KBt^ro)^ - P(ro < t). 

where tq is the hitting time of zero for the Brownian motion B. 

Let ij/^ denote the local time at 1/i^ for the process B./^r^- If we take 
derivatives with respect to t in the equation above, we get 



h'{t) = jE{l - KBt^^,)+ - jP{t, < t) 
K d „ji/K d 



(17) 



The second equality above is due to the Tanaka formula. 

Now to compute the first term on the right of (jl7p . we use [33] Ex. 1.12, 
p. 407] and the explicit (see [2'6\ page 97]) transition function of Brownian 
motion absorbed at zero to get 



g-(l-l/X)2/2t _ ^-{l+l/Kf/2t 
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The second term on the right side of (jl7p above is the density of the first 
hitting time of zero, which we know ([231 pa-ge 80]) to be (27rf^)~-'^/^e~^/^*. 
Thus, we get 



h'it) 



K 



-{l-l/Kf/2t _ -{l+l/KY/2t 



2V27rt 

Thus, h'{t) < if and only of 
2 



1 



-l/2t 



(18) 



Kt 



> e 



l/2t 



^-{\-\IKfl2t _ g-(l+l//^)V2t 





r(2K-i)i 




(2if + l)" 


exp 


2KH 


— exp 


2KH 



(19) 



If i^T > 1/2, it is clear that the above inequality cannot be true in a neigh- 
borhood of zero. Thus, the function h{t) must be initially increasing. 

To get an estimate of the point when it starts decreasing, let t = 1/y. 
Consider the function on the right side of the last inequality. We consider 
two separate cases. First suppose K > 1/2. Then both 2K — 1 and 2K + 1 
are positive. If for two positive parameters A2 > Ai > 0, we define a function 
q by q{y) = exp(Aiy) — exp(— A2y), y > 0, it then follows that 



q'{y) = Xie^'y + X2e-'^'y, 



/(O) = Ai + A2 



q"{y) = Xy^y - Xle-^^y. 
Note that q"{y) < 0, for all 

21og(A2/Ai) 



0<y< 



Ai + A2 



(20) 



(21) 



Since q'{y) is always positive, it follows that q is an increasing concave 
function starting from zero in the interval given by (j2ip . Thus, 



q{y) = q{y) - 9(0) < yq\0) = y (Ai + A2) . 



(22) 



Take Ai = {2K-1)/2K^ and A2 = {2K + 1)/2K'^. Then A1 + A2 = 2/K. 



By dH]) we get that if 



y<Ci ■.= Klog 



then, from (j22|) it follows 





l{2K-l)yl 




jexp 


2K2 


— exp 



2K + 1 
2K -V 



i2K + l)y 
2K^ 



< 2y. 
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That is, by p^ . h'{t) < 0, i.e., h is strictly decreasing for all 



The case when < K < 1/2 can handled similarly. Suppose < Ai < A2 
are positive constants. Consider the function 

r(y) = -Aiy + Asy - e'^iJ' + e'^^^ y G [0, 00). 

Then r(0) = 0, and 

r'{y) = -Ai (1 - e-^^^) + A2 (l - e"^^^) > 0, ye [0, 00), 

because Ai < A2. Thus, for all positive y, we have r{y) > 0, i.e., 

g-Aij, _ g-A2j/ ^ ^_x^+X2)y. 

We use this for Ai = (1 - 2K)/2K^ and A2 = (1 + 2K)/2K^. Note that, as 
before (-Ai + A2)y = 2y/K. 

From (jl9p it follows that > for all t € (0, 00). Thus we have 

established that if K < 1/2, the function t i-> E{l/Xt — K)'^ is strictly 
decreasing for all t G (0, 00). This completes the proof of the proposition. □ 

For mathematical completeness we show below that a similar result can 
be proved for the Bessel process starting from zero, although in this case 
there is no dependence on K. The proof is much simpler and essentially 
follows by a scaling argument. Note that, even in this case the reciprocal of 
the Bessel process is well-defined for all times except at time zero. Hence 
1/Xt, t £ (0,00), can be thought as a Markov process with an entrance 
distribution, i.e., a pair consisting of a time- homogenous Markov transition 
kernel {Pt}, t > 0, and a family of probability measures {/U^}, s > 0, satis- 
fying the constraint jig* Pt = Pt+s ■ Here * refers to the action of the kernel 
on the measure. 

Proposition 9. Let Xf be a 3 -dimensional Bessel process, BES(3), such 
that Xq = 0. For any two time points u > t > 0, and for K > 0, one has 
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Proof. Fix u > t. Recall that BES(3), being the norm of a three dimensional 
Brownian motion, has the Brownian scaling property when starting from 
zero. That is to say, for any c > 0, 

■^Xcs, s>0 j = (X„ s>0), 

where the above equality is equality in law. 

Take c = u/t, and apply the above equality for Xg when s = t, to infer 
that c~^/'^Xu has the same law as Xt, and thus 

Note that for any a > 1, we have (x — aK)^ /a < (x — K)^ , V x > 0. 
Since c > 1, taking a = -y/c, one deduces from ([2] 



e( —-k\ <e(— -K 

\Xu J \Xt 
which proves the result. □ 

We conclude this subsection with an example of a strict local martingale 
S where E(St — K)~^ is not asymptotically decreasing for any K. This, 
coupled with the earlier Bessel result, establishes the fact that functions 
which are not sublinear at infinity can display a variety of characteristics 
when applied to strict local martingales. 

We inductively construct a process in successive intervals + 1) by the 
following recipe. The process starts at zero. The process in the odd interval 
[2i, 2i + l) is an exponential Brownian motion exp(Bt — t/2) starting from 5*21 
and independent of the past. On the even intervals [2i + l, 2i + 2) the process 
5 is an inverse Bessel process starting from S2i+i and again independent of 
the past. The constructed process is always a positive local martingale. The 
value of the function E(St — K)~^ is increasing in the odd intervals due to 
the martingale component, and decreasing (at least when K < 1/2) in the 
Bessel component by Proposition [8j 

One might object to the fact that this process is not strict local through- 
out. But, one can mix the two components, by a sequence of coin tosses 
which decides whether to use Brownian or the Bessel component in the cor- 
responding interval. By choosing the probability of heads in these coins in 
a suitably predictable manner, we can generate a local martingale which is 
strict throughout but E(St — K)~^ does not decrease anywhere. 
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5.1 A multidimensional analogue by Kelvin transform 

In the last section we saw that for any strict local martingale law R there is 
a true martingale law Q such that equality (jl2p holds. The transformation 
g{x) = xh{\/x) has a well-known analogue in dimensions higher than two 
called the Kelvin transform. In our final subsection we present an interesting 
multidimensional generalization of our results for dimensions d > 2. We take 
the following definition from the excellent book on harmonic function theory 
[U Chapter 4] by Axler, Bourdon, and Ramey. 

The Kelvin transform K is an operator acting on the space of real func- 
tions u on a subset of M'^VjO}. Let n be a function on an open subset P 
of M'^\{0}. Let T?* be the image of P under the inversion map 

X ^ X* = x/ . (25) 

For such a u, we define a function K[u] : T>* — t- M by the formula 

K[u]iy) = \y\'-\[y/\y\^). (26) 

Notice that K is its own inverse. 

The most striking property of this transform is that K commutes with 
the Laplacian ([1, page 62]). Let v be the function v{x) = \x\^ Au{x), x G P. 
Then, at any point y G 2?*, we have 

AK[u]iy) = K[v]iy). (27) 

In particular, if u is harmonic in D (i.e., Au = 0), then K[u] is harmonic in 
T>* . Also, if u is subharmonic (i.e.. An > 0), then so is i^[u]. 

A d-dimensional conformal local martingale is a process {X{1), . . . ,X(d)) 
such that each coordinate X{i) is a local martingale and 

{X{i),X{j)) = (X(l)) for ah 1 < i,j < d. 

We consider it to be strict local if at least one of its coordinate processes is 
a strict local martingale. We replace the condition of nonnegativity of one 
dimensional processes by restricting our multidimensional processes in the 
complement of a compact neighborhood of zero. 

Proposition 10. Let D be a compact neighborhood of the origin and denote 
its complement M.'^\D by D'^. Let P be a probability measure on the canonical 
sample space C"^[0,oo) such that, under P, the coordinate process X is a 
conformal local martingale in (Xq = xq € D'^) absorbed upon hitting the 
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boundary of D. Then there is a probability measure Q such that, under Q, 
the coordinate process is a true conformal martingale which takes values in 
(D^y (see (|25p ) such that, for any bounded measurable function U : D'^ — )• 
M, and for any bounded stopping time r, we have 



Fol 



\Xr\^''^U 



Xr 

\XA' 



(28) 



The proof requires the following lemmas. 



Lemma 11. Let |-| denote the Euclidean norm in dimension d. Let ti be 
the hitting time of D i.e., 

Tl = inf {t > : Xt£D] . 
Then, the process \Xt/\T]_\^~'^ , i > 0, is a P-martingale. 

Proof. The function jxl^"*^ is harmonic in M*^. Thus jXfl^^'^ is a local mar- 
tingale itself. Since it is bounded in D^, it must be a true martingale. □ 

Suppose Xq = xq ^ . We can change the law of X by using 1^"*^ 
as a Radon-Nikodym derivative (after normalizing). We get the following 
lemma. 

Lemma 12. Suppose Xq = xq such that xq € D'^. We change P by using 
the positive martingale 

as a Radon-Nikodym derivative. Call this measure Q* . Under Q* , the pro- 
cess 

Y, = (29) 

\XtATi I 

is again a d- dimensional conformal local martingale such that every coordi- 
nate process Y is a true martingale. 

Proof of Lemma\12l Note that, since D is a neighborhood of zero, the set 
[D'^Y is compact. Thus, Y lives in a bounded set. To show that the ith 
coordinate process Y[i) is a local martingale, we use the harmonic function 

u{x) = x{i) on Vi = M"\{0}. 

By (|27p . its Kelvin transform is also harmonic. Hence the process 

\Xt\'-''u{Xt/\Xt\') = \Xt\-'' Xt{i) 
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is a local martingale under P. But, since (until ri) is used as a 

Radon-Nikodym derivative after being scaled, by Bayes rule [231 page 193], 

the process u (^Xt/\Ti/ \Xt/\Ti\'^^ is a local martingale under the changed mea- 
sure Q* . But this implies that Y{i) is a local martingale under Q* . But, 
since every Y{i) is bounded they must be true martingales. 

To show that Y is conformal, we use the harmonic function u{x) = 
x{i)x{j) again on the full domain M"\{0} for any pair of coordinates i 7^ j. 
Exactly as in the previous paragraph, we infer from ()27p that Yt{i)Yt[j) is 
a local martingale under Q. But this implies (Y{i),Y{j)) = 0. That their 
quadratic variations must be the same follows from symmetry. □ 

Proof of Proposition [TOl The construction of Q has been done in Lemma 
[T2] where it is the law of the process Y under Q* . Note that the Radon- 
Nikodym derivative (i){Xt) never hits zero under P. Thus, even under Q* , 
the process Y never hits zero. To show the equality (j28|) . we use the change 
of measure to get 

Xr\^~'^ \Xr\'^~'^ U {Xr) 
= \x^\'''^EP [U{Xr)]. 

This completes the proof of the result. □ 

6 Applications to financial bubbles 

A natural question is: what happens to a financial market when the no 
arbitrage condition yields a strict local martingale (rather than a true mar- 
tingale) under a risk neutral measure? Several authors have looked at this 
problem and offered solutions to anomalies which might result from the lack 
of the martingale property. One interesting perspective offered in this di- 
rection is the theory of price bubbles as argued in 2000 by Loewenstein and 
Willard l25j. They propose that to identify a bubble one needs to look 
at the difference between the market price of an asset and its fundamental 
price. Their argument is later complemented and further developed by Cox 
and Hobson [4] and the two articles by Jarrow, Protter, and Shimbo |20j . 
|21j . In particular they identify that a stock price which behaves as a strict 
local martingale under an equivalent local martingale measure is an example 
of a price bubble. Moreover, Cox and Hobson [4J go on to exhibit (among 
other things) how in the presence of bubbles put-call parity might not hold 
and call prices do not tend to zero as strike tends to infinity. 



\YJ'^~'^U 



Yr 
\Yr? 



I \d-2 TT^P 

|xo| E 
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We consider a market with a single risky asset (stock) and zero spot 
interest rate. Let {St}, t G (0,c«), be a positive continuous strict local 
martingale which models the discounted price of the (non-dividend paying) 
stock under an equivalent local martingale measure. We have the following 
result which follows immediately from Corollary [T] and the subsequent Bessel 
example. 

Proposition 13. Suppose for a European option, the discounted pay-off at 
time T is given by a convex function /i(S'r) which is sublinear at infinity, 
i.e., lim^^^oo h{x)/x = 0. Then the price of the option is increasing with the 
time to maturity, T , whether or not a bubble is present in the market. In 
other words, E{h{ST)) is an increasing function ofT. For example, consider 
the put option with a pay-off {K — x)'^ . 

However, for a European call option, the price of the option E[St — K)'^ 
with strike K might decrease as the maturity increases. 

This feature may seem strange at first glance, but if we assume the 
existence of a financial bubble, the intuition is that it is advantageous to 
purchase a call with a short expiration time, since at the beginning of a 
bubble prices rise, sometimes dramatically. However in the long run it is 
disadvantageous to have a call, increasingly so as time increases, since the 
likelihood of a crash in the bubble taking place increases with time. 

It is worthwhile to compare our results with processes that are increasing 
in the convex order. We refer the interested reader to the articles by Hirsch 
and Yor |17j . |18j and the references within. 

Of course, pricing a European option by the usual formula when the 
underlying asset price is a strict local martingale is itself controversial. For 
example, Heston, Loewenstein, and Willard [16] observe that under the ex- 
istence of bubbles in the underlying price process, put-call parity might not 
hold, American calls have no optimal exercise policy, and lookback calls have 
infinite value. Madan and Yor [27] have recently argued that when the un- 
derlying price process is a strict local martingale, the price of a European call 
option with strike rate K should be modified as lim^^oo E {{Stat^ ~ K)^], 
where r„ = inf {t > : > n}, n G N, is a sequence of hitting times. This 
proposal does however, in effect, try to hide the presence of a bubble and act 
as if the price process is a true martingale under the risk neutral measure, 
rather than a strict local martingale. 

A different approach to such market anomalies has been studied exten- 
sively by several authors when the candidate Radon-Nikodym derivative for 
the risk-neutral measure turns out to be a strict local martingale. This is the 
spirit taken in Platen [31], Platen and Heath |32j, Fernholz and Karatzas 
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|12j . Fernholz, Karatzas, and Kardaras |13j . and Karatzas and Kardaras 
|22j . For the latter group of authors this is intimately connected with the 
concept of a weakly diverse market and displays a number of anomalies sim- 
ilar to the case of bubbles. For example, put-call parity fails to hold in such 
markets. See, Remark 9.1 and 9.3, and Example 9.2 in [13j. Finally, we 
mention that P. Guasoni and M. Rasonyi have recently studied the robust- 
ness of bubble models by showing that if is a price process which models 
a financial bubble, under some conditions one can find another process 
which is "e close to S" and is not itself a bubble. Therefore in this sense 
models for bubbles perhaps can be seen as delicate from the standpoint of 
model risk [15]. 
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